A fluctuation theorem for phase turbulence of chemical oscillatory waves 
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Through numerical simulations of the Kuramoto equation, which displays high-dimensional dis- 
sipative chaos, we find a quantity representing the cost for maintenance of a spatially non-uniform 
structure that appears in the phase turbulence of chemical oscillatory waves. We call this quan- 
tity the generalized entropy production and demonstrate that its distribution function possesses a 
symmetry expressed by a fluctuation theorem. We also report a numerical result which suggests a 
relation between this generalized entropy production rate and the Kolmogorov-Sinai entropy. 

PACS numbers: 05.45.+b,82.40.Bj 



When a system exists in a non-equilibrium steady 
state, the entropy produced in the system continuously 
flows into the environment. The amount of the entropy 
produced during a finite time interval is called the 'en- 
tropy production'. As stipulated by the second law of 
thermodynamics, the entropy production must be posi- 
tive. 

Entropy production can be interpreted as the cost for 
the maintenance of non-equilibrium steady states. This 
leads us to expect that we can find a quantity which can 
be regarded as a generalized entropy production charac- 
terizing non-thermodynamic as well as thermodynamic 
systems. With the hope of demonstrating this point, we 
study the phase turbulence of chemical oscillatory waves 

i- 

When a pacemaker is situated in a small region in a 
spatially uniform phase turbulent state, a spatially non- 
uniform state appears, as is easily understood. In this 
paper, through numerical experiments, we find an ex- 
pression for a generalized entropy production character- 
izing such a spatially non-uniform state. We further at- 
tempt to relate this generalized entropy production with 
the Kolmogorov-Sinai (KS) entropy, which measures the 
rate of information loss in chaotic dynamical systems. 

We first introduce a mathematical model describing 
phase turbulence. For simplicity, we consider reaction 
diffusion systems in a one-dimensional circuit. Such sys- 
tems are considered as describing the behavior of the con- 
centrations of chemical species. For a system of this type, 
when certain conditions are satisfied (see § for details), 
their concentrations oscillate at each position x, and the 
phase <f> of the oscillation varies slowly in time t. In the 
weakly unstable case, the time evolution of the phase 
qb{x, t) is described by the Kuramoto equation |jl|, 0], 



d<t> d 2 - 
dt dx 2 



dx 4 



d<p_ 



(i) 



where we have carried out appropriate scalings for x, t 
and 4> so as to obtain the simplest form of the equation. 
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The Kuramoto equation (|lj) is defined over the region 
< x < L, with periodic boundary conditions. The 
right-hand side of ([[]) represents the effect of external 
oscillation introduced by the pacemaker: e corresponds 
to the pacemaker strength, and the form of uj{x) reflects 
its position. We assume that lu(x) = 1 in the region 
< x < H and oj(x) — otherwise. In the analysis 
below, we study the case (L, H) = (128, 10). 

We obtain an expression for the generalized entropy 
production by undertaking a numerical analysis of (lj). 
In our numerical simulations, we discretized space with a 
unit size Sx = 1.0 and employed an explicit Euler method 
with time step St = 0.01. The nonlinear term was dis- 
cretized as ((pi+i ~ Pi~i)/25x) 2 /2, where pi — <j>(i6x,t). 
We believe that this numerical scheme is sufficient to al- 
low for investigation of the behavior of interest. 

When there is no pacemaker (e = 0), (|]) possesses 
spatio-temporal chaotic solutions for a rather wide class 
of initial conditions j3, H|. (Precisely speaking, there is 
a very small basin for stable spatially periodic solutions 
|Q, ||, 0|. However, we do not consider such solutions, 
because this basin is too small to be observed when ini- 
tial conditions are assigned randomly 0.) The statistical 
properties of this spatio-temporal chaos have been exten- 
sively studied both numerically j^, [|, Q and theoretically 
II- 

When e ^ 0, the system possesses spatially non- 
uniform statistical properties. As an example, in Fig. 
[I], we plot the long-time average of the phase profiles, 
each of which is shifted so as to satisfy cf>(0, t) = 0. 

Now, let us consider the generalized entropy produc- 
tion which represents the cost for maintaining the spa- 
tially non-uniform state. It is important to note here that 
the generalized entropy production will not be expressed 
in terms of heat, in contrast to the usual entropy produc- 
tion of non-equilibrium steady states, because the phase 
turbulent state cannot be described by quantities char- 
acterizing thermodynamic equilibrium states. For this 
reason, we have no intention to base our analysis of the 
generalized entropy production on thermodynamic con- 
siderations. Instead, we consider a more abstract condi- 
tion that we expect the generalized entropy production 
to satisfy. 
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FIG. 1: Long time average of the phase profiles, each of 
which is shifted so as to satisfy (j>(0,t) = 0. Here, e — 0.1. 
The averaging time was 4 x 10 7 time units for the statistical 
steady state. 



FIG. 2: logn(z;r, 1) -log n(-z; r, 1). The dotted line corre- 
sponds to log II(z; r, l)-logII(— z; r, 1) = 0.242. Here, e = 0.1 
and r = 200. The distribution was made by using 2 x 10 5 tra- 
jectory segments for the statistical steady state. 



As a recent development in the study of non- 
equilibrium steady states 0, [□], [H| an interesting 
relationship regarding the entropy production has been 
established. Defining H(z;t) to be the probability den- 
sity that the entropy production S T during a finite time 
interval r takes the value z in a given non-equilibrium 
steady state, it has been shown that this probability den- 
sity possesses the symmetry 



log Il(z; r) — log II(— z; r) = z + o(r) 



(2) 



in the limit r — > oo. This symmetry is referred to as a 
'fluctuation theorem'. From (^), it is easily derived that 
the most probable value of the entropy production ratio 
S r /r in the large r limit is positive. Considering the 
above result for entropy production in non-equilibrium 
steady states, we conjecture that a similar result can be 
obtained for a wider class of nonequilibrium states by 
appropriately defining a generalized entropy production 
and seeking a generalized fluctuation theorem it satisfies. 

The problem we face is to find a generalized entropy 
production characterizing spatially non-uniform phase 
turbulent states, whose distribution function satisfies 
(||). By recalling that the entropy production in a non- 
equilibrium steady state can be expressed in the form 
[force] [flux] /[temperature], and by interpreting ed x u and 
d x <p as corresponding to a force and a flux, we regard 
the following quantity as a candidate of the generalized 
entropy production: 



e/3 



dt 



L Q2 
dx( — Ld(x)) 



dx 2 ' 



(3) 



where we note that the integrand in is rewritten as 
d x u>d x 4>. corresponds to the effective inverse tempera- 
ture of the phase turbulence whose value is determined 
later so that the distribution function of £55 satisfies (0) . 
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FIG. 3: Entropy production rate a as a function of e. The 
dotted line corresponds to a — e 2 . The averaging time was 
4 x 10 time units for the statistical steady states. 

Through numerical experiments, we obtained the dis- 
tribution function H(z; r, (3) of for a statistically 
steady state. Since the value of f3 is not determined 
yet, we first considered an arbitrary value and calculated 
II(z; r, 1). As seen from Fig. |^, it seems that the equality 

logn(z;r, 1) -logU(-Z]T, 1) = 0.24z (4) 

holds approximately. Then, since the trivial identity 

H(z;T,0)=IL{z/piT,i)/0 (5) 

holds, r, (3) with = 0.24 satisfies the fluctuation 
theorem (||). 

If we wish to characterize spatially non-uniform phase 
turbulence by a single quantity, the generalized entropy 



3 



production ratio, 



a = lim — — , 

T — >oo T 



(6) 



may be most important. In Fig. [|, we plot a as a function 
of the pacemaker strength e. As is true generally for non- 
equilibrium steady states, we see that a is proportional 
to e 2 in the region of small e. 

We can formulate an argument to support our nu- 
merical result by considering the Yakhot conjecture 
H [TtJ [l^l that the long time and large distance be- 
havior of statistical properties of chaotic solutions to (|l|) 
are equivalent to that of solutions to the stochastic evo- 
lution equation (which is essentially the same as the Bur- 
gurs equation with noise) 



~\ 2 




(7) 



where £ is a Gaussian white noise satisfying 

(£(x, t)£(x', t')) = 2DS(t - t')S(x ~ x'). (8) 

Here, v and D are positive constants. It is quite easy 
to demonstrate a fluctuation theorem for this stochastic 
system (See Refs. ^3|). The result is the following: 
When we define the entropy production as 



D 



dt 



L q2 

dx(—uj(x)) 



dx 2 1 



(9) 



the distribution function of for a statistical steady 
state satisfies the fluctuation theorem (||). Therefore, 
if the Yakhot conjecture is valid and if this conjecture 
ensures the equality of the space-time inetegration in (|^) 
with that in (^), the expression of (||) should hold with 
(3=v/D. 

We would like to consider the generalized entropy pro- 
duction from the dynamical system viewpoint. Until 
now, we have not succeeded in deriving the fluctuation 
theorem in a manner based on dynamical system consid- 
erations. However, we have obtained a numerical result 
concerning the relation between the generalized entropy 
production rate and the KS entropy, which can be calcu- 
lated using the Lyapunov analysis. We now discuss this 
result. 

Performing numerical integrations of the linearized 
equation corresponding to (]l|) for chaotic solutions and 
employing the Gram-Schmidt decomposition technique 
[jlif , we calculated the long time average of the expansion 
ratio, A(m), for m-dimensional volume elements along 
chaotic trajectories 0. In Fig. 4, the graph of A(m) is 
plotted for the system without the pacemaker. The KS 
entropy is defined as the maximal value of A(m). Let 
h^{e) be the KS entropy for the system with pacemaker 
strength e. Since we naively expect that <r(e) is related 
to the quantity 



A/i ks (e) = h ks (e) - h ks (0), 



(10) 




FIG. 4: Long time average of the expansion ratio of m- 
dimensional volume elements. The averaging time was 10 4 
time units for the statistical steady state. 




FIG. 5: (a, -A/i ks ) for different values of e (0 < e < 0.2). 
The dotted line corresponds to — Ah^ B — a/2. The averaging 
time was 4 x 10 time units for the statistical steady states. 



we plot points (a, 
displayed in Fig. F> 



-A/iks) for different values of e. 
the results indicate the relation 



As 



A/l ks = \(T 



(11) 



for small e. Note that we observed a greater deviation 
from (^) when e is increased. 

We do not understand why the numerical factor 1/2 
appears in (]lT|). However, the fact we have obtained such 
a simple relationship for the system we consider presently 
may suggest the existence of a universal relation between 
the generalized entropy production and some dynami- 
cal system quantity. (In the study of a time-dependent 
Hamiltonian system, the same factor appears in the rela- 
tion between the Boltzmann entropy difference and the 
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excess information loss found through Lyapunov analy- 
sis Further investigation is necessary to clarify this 
point. 

We now make some remarks. 

A fluctuation theorem was first demonstrated numer- 
ically [jl0| and proved mathematically jll| in the study 
of a deterministic Gaussian thermostatted model. In this 
model, the entropy production is given by the phase space 
contraction, and the fluctuation theorem is derived by 
exploiting the time-reversibility of the deterministic evo- 
lution equation. Indeed in many cases, including that of 
the stochastic model ^ with (||), the entropy produc- 
tion for a given trajectory has been related to the ratio 
of the measure of this trajectory to that of its time re- 
versal trajectory. However, such an argument cannot be 
applied to ([!]), because it does not possess time-reversal 
symmetry. 

On a more abstract level, as discussed by Maes [ p"3| , 
fluctuation theorems can be associated with a transfor- 
mation P, that acts on trajectories and satisfies P 2 = 1. 
In certain situations (e.g. for steady states), this trans- 
formation is equivalent to time reversal, but in general 
this is not the case. The problem is thus to determine the 
form that P takes in the present case and to determine 
how to derive our fluctuation theorem from it. 

Finally, we briefly discuss the size dependence for the 
system we have studied. It is known that the statisti- 
cal properties of phase turbulence without a pacemaker 
depend on L in an anomalous way (|, |9[ [l6| [l7[ [l8| . 
The system size in our numerical experiments is much 



smaller than the size representing cross-over to the region 
in which the anomalous behavior is observed. It is an in- 
teresting problem to extend our study to this anomalous 
scaling region. 

In summary, we have demonstrated that the gener- 
alized entropy production (^|), which has been obtained 
phenomenologically in the consideration of spatially non- 
uniform phase turbulent states, satisfies the fluctuation 
theorem (^). Unfortunately, there is no obvious way to 
derive this fluctuation theorem from the deterministic 
evolution equation in the present case, because we cannot 
employ an argument relying on time-reversibility. Nev- 
ertheless, using the Yakhot conjecture, we have formu- 
lated a reasonable explanation of the result. We have 
also obtained the numerical result (fill), which suggests 
an interesting relation between the generalized entropy 
production ratio and the KS entropy. We believe that 
a deeper understanding of the problem considered here 
will allow for the study of dissipative high-dimensional 
systems, such as fluid turbulence and granular flow, from 
a new point of view. 
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